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REGULARIZATION STRATEGY FOR INVERSE 
PROBLEM FOR 1+1 DIMENSIONAL WAVE 
EQUATION 

JUSSI KORPELA, MATTI LASSAS AND LAURI OKSANEN 


Abstract. An inverse boundary value problem for a 1+1 dimen¬ 
sional wave equation with wave speed c(x) is considered. We give 
a regularisation strategy for inverting the map A : c i—A, where A 
is the hyperbolic Neumann-to-Dirichlet map corresponding to the 
wave speed c. More precisely, we consider the case when we are 
given a perturbation of the Neumann-to-Dirichlet map A = A + £", 
where £ corresponds to the measurement errors, and reconstruct 
an approximate wave speed c. We emphasize that A may not not 
be in the range of the map A. We show that the reconstructed 
wave speed c satisfies ||c — c||loo < C||if||^A8^ regularization 
strategy is based on a new formula to compute c from A. 
Keywords: Inverse problem, regularization theory, wave equa¬ 
tion. 


1. Introduction 

We consider an inverse bonndary valne problem for the wave eqnation 
(dt -c(xfd^)u(t,x) = 0. 

We introdnce a regnlarization strategy to recover the sonnd speed c(x) 
by nsing the knowledge of pertnrbed Nenmann-to-Dirichlet map A. 
Onr approach is based on the Bonndary Control method [21 |6l l53] . 

A variant of the Bonndary Control method, called the iterative time- 
reversal control method, was introdnced in |2]. The method was later 
modified in na to focns the energy of a wave at a fixed time and in jlB] 
to solve an inverse obstacle problem for the wave eqnation. Here we 
introdnce yet another modification of the iterative time-reversal control 
method that is tailored for the 1+1 dimensional wave eqnation. 

Classical regnlarization theory is explained in [16]. Iterative regnlar¬ 
ization of both linear and nonlinear inverse problems and convergence 
rates are discnssed in Hilbert space setting in [lOl UHl [20l HU 03] and 
in Banach space setting in [I9l |23l |2ll [301 HU HU HSj. Onr new resnlts 
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give a direct regularization method for the nonlinear inverse problem 
for the wave equation. The result contains an explicit (but not neces¬ 
sarily optimal) convergence rate. 

By direct methods for non-linear problems we mean explicit con¬ 
struction of a non-linear map solving the problem without resorting 
to a local optimisation method. In our case the map is given by (55). 
The advantage of direct approaches is that they do not suffer from the 
possibility that the algorithm converges to a local minimum. In par¬ 
ticular, they do not require a priori knowledge that the solution is in 
a small neighbourhood of a given function. There are currently only 
few regularized direct methods for non-linear inverse problems. An 
example is a regularisation algorithm for the inverse problem for the 
conductivity equation in |31| . Also, a direct regularized inversion for 
blind deconvolution is presented in 


1.1. Statement of the results. Let Co,Ci, L,m > 0 and dehne the 
space of velocity functions 

(1) V(A) :={c e L°°{M)-Co < c{x) < Ci, 

II'"IIc2(m) — ^ ~ ^ ^ ^0 (0) 

where we denote by M the half axis [0, cxd). Let 

(2) T>T. 

Uo 

For c G T>{A) and / G L^(0,2T), the boundary value problem 

(3) {d^ — c{x)'^dl)u{t,x) = 0 in (0,2T) x M, 
d^u{t,0) = /(f), 

u\t=o = dtu\t=o = 0 , 

has a unique solution G H^{{0,2T) x M). Using this solution we 
dehne the Neumann-to-Dirichlet operator 

(4) A : L\0,2T) ^ L\0,2T), Af = u^,=o- 
We dehne for a Banach space E 

C{E) := {A : E ^ E] A is linear and continuous}. 

Let X = L°°{M) and Y = £(L^(0,2T)). We dehne the direct map 

(5) A : V{A) CX A{c) = A. 

We show in Appendix A, Theorem that the maps (|^ and (|^ are 
continuous. Here the range Ran{A) = A{T>{A)) and the domain ^{A) 
are equipped with the topologies of Y and X, respectively. 
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We consider the inverse problem to recover the velocity fnnction c 
by using the boundary measurements A. It is well-known that A is 
invertible. Let us record the following: 

Theorem 1. The inverse map 

A-^ : Ran{A) C F -> V{A) C X, = c, 


is continuous. 

For the convenience of the reader we give a proof of Theorem in 
Section 2, where we also give a new formula to compute c from A. 
Our main result concerns perturbations of the Neumann-to-Dirichlet 
operator of the form 


( 6 ) 


A = A + f. 


where £ & Y models the measurement error. We assume that ||T||y < 
e, where e > 0 is known. In this situation we can not use the map 
A~^ to calculate function c since A may not be in the range TZ{A). We 
recall the definition of a regularization strategy, see e.g. [12] and |3Uj . 

Definition 1. Let X,Y he Banach spaces and ^{A) C X. Let A : 
BiA) Y be continuous mapping. Let ao G (0, oo]. A family of 
continuous maps TZa : F —)■ X parametrized by 0 < a < ao is called a 
regularization strategy if 


lim 7la(A(c)) = c 

a^O 


for every c G V{A). A regularization strategy is called admissible, if 
the parameter a is chosen as a function of the noise level e so that 
IhUe^o Q:(e) =0 and for every c G ViA) 


lim sup 

e-5-O 


'R'a{e)A. C 


: A G F, 

X 


A — A{c) 



0. 


Figure gives ns a schematic illustration of regularization. 
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Model space Data space 



Figure 1. The idea of regularization is to construct a 
family TZa(e) of continuous maps from the data space Y 
to the model space X in such a way that c can be approx¬ 
imately recovered from noisy data A. For a smaller noise 
level e the approximation 7la(e){A.) is closer to c. More 
details and a similar figure in a more general setting can 
be found in |1H Fig. 11.5]. 


The main result is presented next and says that we have an admis¬ 
sible regularization strategy inverting A. 

Theorem 2. There exists an admissible regularization strategy TZ^ with 
the choice of parameter 

/ N 13„4 4 

a{e) = 2 9 Tses 


that satisfies the following: For every c G T>{A) there is eo, C > 0 such 
that 


sup 


'hZa{t)^ C 


X 


: A EY, 


A{c) — A 


Y 


< e ^ < Ceis, 


for all e E (0, Cq). 


We will give explicit choices of TZa and cq, see (54) and (55) below. 


1.2. Previous literature. From the point of view of uniqueness ques¬ 
tions, the inverse problem for the 1+1 dimensional wave equation is 
equivalent with the one dimensional inverse boundary spectral problem. 
The latter problem was thoroughly studied in 1950s dziEaiig and we 
refer to [221 pp. 65-67] for a historical overview. In 1960s Blagovescen- 
skii dans] developed an approach to solve the inverse problem for the 
1+1 dimensional wave equation without reducing the problem to the 
inverse boundary spectral problem. This and later dynamical methods 
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have the advantage over spectral methods that they require data only 
on a finite time interval. 

The method in the present paper is a variant of the Boundary Control 
method that was pioneered by M. Belishev |2] and developed by M. 
Belishev and Y. Kurylev mu in late 80s and early 90s. Of crucial 
importance for the method was the result of D. Tataru |53j concerning 
a Holmgren-type uniqueness theorem for non-analytic coefficients. The 
Boundary Control method for multidimensional inverse problems has 
been summarized in ElES], and considered for 1+1 dimensional scalar 
problems in lac! and for multidimensional scalar problems in |251 l28l 
l33l EHl EH] For systems it has been considered in |3ll ESj. Stability 
results for the method have been considered in jl] and [29]. 

The inverse problem for the wave equation can be solved also by using 
complex geometrical optics solutions. These solutions were developed 
in the context of elliptic inverse boundary value problems [52], and in 
|45j they were employed to solve an inverse boundary spectral problem. 
Local stability results can be proven using (real) geometrical optics 
solutions PEHIEI], and in BOI a local stability result was proved by 
using ideas from the Boundary Control method together with complex 
geometrical optics solutions. Finally we mention the important method 
based on Carleman estimates na that can be used to show stability 
results when the initial data for the wave equation is non-vanishing. 


2. Modification of the iterative time-reversal control 

METHOD 


In this section we prove Theorem in such a way that we can utilize 
the proof to construct a regularization stratezy as in Theorem Let 
A be as defined in (|^. Let r G [0,T]. We define linear operators in Y 
by 

(7) 7A(^,s)/(s)ds, 

Rf{t) = f{2T -t), K = JA- RARJ, 

Bf{t) = 1(0,T)(^) f{s)ds, Prf{t) = l(o,r)(t)/(t), 


where 


1a (A s) 


1, t + s < 2T and s > f > 0, 
0, otherwise 
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and 


1(0,r) (t) 


1, te(0,r), 

0, otherwise. 


Let / G L^(0,2T). Using solntion G if^((0,2T) x M) of Q we 
define 

(8) Ut ■■ L‘^{0,2T) ^ H\M), Ur/ = nl(T). 

Let ns denote dV = c~‘^dx and recall the Blagovestchenskii identities 

(9) {u'^{T)A)L^{M-,dv) = (/,-Bl)L2(o,2r), 

(m1(T), = {f ■,Kh)i2(^Q2T)- 

The identities (|^ originate from the work by Blagovestchenskii [11] 
and their proofs can be fonnd e.g. in jU]. We dehne the domain of 
inflnence 


( 10 ) 


M{r) = {x G M; d(x, 0) < r}, 


where d{x, 0) = We nse the following resnlt that is closely 

related to [9] 

Theorem 3. Let r G [0,T] and a > 0. Let K,B, and he as defined 
in &■ Let us define 

(11) Sr = {fe L\0,2T) : snpp(/) C [T-r,T]}. 

Then the regularized minimization problem 

(^(/) -^/)l2(o,2T) ~ 2{f, B1 )l2(q 2T) + « II/IIl2(o,2T)) ’ 
has unique minimizer 

(12) fa,r = {PrKPr + a)~^ PrBl 

and the map r i—)■ fa,r is continuous [0,T] —)• L^(0,2T). Moreover 
uP’^iT) converges to the indicator function of the domain of influence, 


mm 

f&Sr 


(13) 


limJ|nl“-(T)-lM(.)l|^.(^.^^) 


= 0 . 


For the convenience of the reader we give a proof. 

Proof of Theorem Let a > 0 and let f E Sr- We dehne the energy 
fnnction 


(14) E{f) := (/, Kf)L^o,2T) - 2{f, B1 )l2^o,2T) + « ||/|| L2(o,2T) • 


REGULARIZATION FOR WAVE EQUATION 


7 


The finite speed of wave propagation implies supp(n'^(T)) C M{r). 

Using (|^ we can write 

(15) 

Ilf l|2 II ||2 2 

f'(/) = Ip {T) — ^M{r)\\ 1^2 ■,(!¥) ~ ln^fi)llL2(M;dV) “ II/II L2 (o,2T)) ' 

Let C Sr be such that 

lim E{fj) = inf E{f). 

j-i-oo /eSj. 

Then 

^ ll/illL2(0,2T)) — + ll^^h) IL2(;v^.(^y) ) 

and we see that is bounded in Sr- As Sr is a Hilbert space, there 

is a subsequence of (/j)^i converging weakly in Sr- Let us denote the 
limit by f^o G Sr and the subsequence still by 

By Theorem 1^ in Appendix A below, the map Ut '■ f •—t {T) is 

bounded 

Ut : L\0,2T) H\M). 

The embedding I : H^{M) ^ L‘^{M) is compact and thus Ut '■ f ^ 
{T) is a compact operator 

Ut : L2(0,2T) ^ L^{M). 

Hence we have a subsequence {fj)^i for which uS(T) —)• u^°°{T) in 
L‘^{M) as j —>■ oo. Moreover, the weak convergence implies 

II/oo|Il2(o,2T)) — lim inf ||/j||i 2 (o, 2 T) • 


Hence 


^^wllL2(jvf.rfy) l|lAf(r)||^2(M.rfy) 


+ Oi ||/( 


oo|Il2(o,2T) 
2 


- ^^h)IL2(M;dy) PA7(r)||^2(^.^y) + «lim^f ll/illz,2(0,2r) 


j^oo 

= liminfL;(/j) = inf E{f), 

J-!-00 f&Sr 


and thus /oo G Sr is a minimizer for (14). We denote by the Frechet 
derivative to direction h. If 


0 — DhE{f) — 2{h, PrKPrf) L2(^Q2T)~‘2{h, PrBl) T2(0^2T)+‘2<^{h, f) L'^(0,2T), 
for all h G L^(0, 2T), then 

{PrKPr + a)f = PrBl. 

Using (|^ we have 

{{PrKPr + a)/,/)n2(0,2T) = (M^''^(T),n^’-^(T))i2(M;,y) + (“/j /)l2(o,2T)- 
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Operator PrKPr + a is coersive when a > 0. The Lax-Milgram The¬ 
orem implies that it is invertible, and we have an expression for mini- 
mizer 


fa,r ■= /oo = {PrKPr "h a) ^PrBl. 


According to j33], see also |27) . we know that 


{u>{T) e D(M(r))i / e Sr] 


is dense in L‘^{M{r)). Let e > 0 and £x G Sr such that 


(16) 


uS{T) - 


2 

L'^{M-,dV) 


< e. 


Using (15) we have 


uK^{T) - lM(r) 


2 

L'^{M-,dV) 


< E{fa,r) + ||l 


M(r) \ 


2 

L^{M;dV) ■ 


Because E{fa^r) < E{f\) we have 


\u 


fa,' 


1 12 Ilf 112 2 

{T) - ^M{r)\\^2i^M-,dV) - Ir “ ^^(p\\L'^{M-,dV) ' 


Using (16) we may choose first small e > 0 and then small a > 0 to see 
that v/°^) tends to Imp) in K{M) as a —?• 0. □ 


See Figure [^for a visualization of M{r). 
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Figure 2. When sending information from an interval 
of length r, that is when supp(/) C [T—r, T], the solution 
u^{x,T) is supported in the domain of influence M{r). 


We dehne the travel time coordinates for x G M hy 
T : [0, cx)) —)■ [0, cx)), t{x) = d{x, 0). 

The function r is strictly increasing and we denote its inverse by 

X = r~^ ■ [0, oo) -)■ [0, cx). 


We have 

(17) X(0) = 0, x'(t) = = c(x(t)). 


Thus denoting v(t) = c(x(t)) and use V(r) to denote the volume of 
M(r) with respect to the measure dV we have 


(18) V(r) 


'•X(r) 


dx 


c(x] 


X'{i)dt 

v{ty 


dt 

v{t)' 


Note that M(r) = [0,x(r)]. In particular, V{r) determines the wave 
speed in the travel time coordinates, 

u(r) = 


(19) 


drV (r) ’ 
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and also in the original coordinates since 


c(x)=v{x ‘(x)), ,\(i) = / 


( 20 ) 


Using Theorem and (|^ we have a method to compnte the volnmes 
of the domains of inflnence 

(21) V(r) = ||lM(r)||i2(M;dU) = 

where r G [0,T]. We are ready to prove Theorem 

Proof of Theorem \1\ For a given measnrement A, Theoremj^and eqna- 
tions (19), (20), (21) give ns a way to calcnlate for all x G (0, L) the 
valne of the velocity fnnction 

c{x) = n(x“^(a:)) = Al"^(A)(a;). 


As we assnmed that ontside of the interval (0, L) the fnnction c is iden¬ 
tically one, the proof for the existence of inverse map A~^ is complete. 
Since the direct map A is continnons and T>{A) is relatively compact 
in X, we see that A~^ is a continnons map. □ 


3. Stability of regularized problem 


In this section we prove Theorem We will constrnct the opera¬ 
tor Tla{£) as a composition of several operators. The constrnction is 
motivated by the proof of Theorem We dehne for a Banach space E 

}C{E) = {A E C{E)] A is compact}. 

Let J, R be as dehned in ([^. Using ([^ we see that J G /C(L^(0,2T)). 
We dehne 

(22) K :Y ^ /C(L2(0, 2T)), KA = JA- RARJ. 

if : r ^ C'([0,T],U), HA = r Pr{KA)Pr. 


Proposition 1. We have ||fi||y_^^QQ < 2T. 

Proof. Let r G [0,T]. We have estimates ||Pr||y < 1, ll-^lly ^ 1; 
ly < T, and 


HA{r) 

Thns 


z:(l2(o,2T)) 


< 2 


£(L2(o,2T)) 


A 


z:(l2(o,2T)) 


< 2T 


A 


£(L2(0,2r)) 


liil 




y[o,r],y) — 


It remains to show that r i—>■ HA{r) is continnons. Let ns denote 
K = KA. Let r, s G [0,T]. We nse the singnlar valne decomposition 
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for the compact operator K. There are orthonormal bases {(j)n}’^=i G 
L^(0,2T) and G L^(0,2T) such that 


(23) 


^ ^ h'n(/) 0n)L^(O,2T)'0n) 


n=l 


for all / G L^(0,2T), where G M are the singular values of K. We 
dehne the family of hnite rank operators by the formula 


(24) 


K'^f — ^/in(/, 0n)L2(O,2r)l/’n- 


n=l 


Then 
(25) ^ 

PrKPrf - P,KPJ 


< 


L2(0,2T) 

PrKPrf - PrK^Prf 


L2(o,2T) 


PsK^Prf - PsK^Psf 


L2(o,2T) 


+ 


PrK'^Prf - PsK^Prf 
PsK^Psf - PsKPsf 


L2(o,2T) 


+ 


L2(o,2T) 


Let e > 0 and let 11/11^2(0 2 r) — choosing m large enough we have 


PrKPJ - PrK^Prf 


L2(o,2T) 


+ 


PsK^Psf - PsKPsf 


e 

< 

L2(o,2T) 2 


Applying projections to (24) we see that 


PsK^Prf = J2f^n{f,Prfn)Psi^r. 


n=l 


For the second term in the sum 


PrK^Prf - PsK^Prf 


n=l 


L2(0,2T) 
m 

— ^ ^ \^^n\ II (-Pr ~ -Ps) 0 n|| L‘^{0^2T) — CijTl) |r — s| 2 . 


we have an estimate 


'^fln{f,Pr(pn){Pr “ Ps)A 


L2(o,2T) 


n=l 


For the third term in the sum we have an analogous estimate 


PsK^Prf - PsK^Psf 


L2(o,2T) 


=< C{m)\r — s| 2 . 


2 

Putting these estimates together and choosing |r — s| < 5(e) = , 

we see that 


PrKPr - PsKPs 


< e. 


□ 


Y 
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Let US define 

(26) Ml = sup{||^(c)||^(^2(0 2r));cer'(^)}. 

Using the continuity of A, see Theorem below, we see that Mi < cx). 
We dehne M 2 = 2TMi. Let c G 'D{A) and denote A = A{c). We use 
again the notations H = HA, H = HA and = HA{r). Using 
Proposition we have 

(27) ll-^llc([o,T],y) — ^ 2 - 

We dehne M 3 = M 2 + 3 and a family }ag(o, 2 ] ^ C(M) by 
( 0 , if s > M 3 - f, 

= l + ‘f s6(a^3-|,Ms-|1, 

[ 1 , if s<M3-f. 

For a G (0,2] we dehne 

(28) Z„:C'([0,T],U)^C'([0,T],U), 

-1 


Z^{H)=r^^i 


M3-(M + a) ]{Hr + a 

c([o,r],r) / V 

Let be a Banach space and let if G -U. Let e > 0. We denote 


(29) 


BE{H,e) ■.= {HeE: 


H-H 


<e}. 


Proposition 2. Let e G (0,1) and let p G (0, |). Let a = 2e^ and 


let \\H\ 


C ([0 


T]Y) — "^ 2 - Let Hr & Y he positive semidefinite. Let us 


assume that H G Bc{[q,t],y){H, e). Then 




c([o,T],y) 


Proof. Using (28) we see that if \1'„( M 3 — {H + a) 
then 

M 3 - (5 + a) 


c([o,r],Y) 


7 ^ 0 , 


a 


< M3 — w < M3 
c([o,r],y) 4 


-1 


and [Hr + a) is dehned by the formula 


-wA 


M 3 - (Hr + a) \ -1 
M 3 


Me 


E 

1=1 


M 3 - [Hr + a) 


M 3 


This gives that Za{H){r) G Y, when r G [0,T]. Proposition 1 


gives 


continuity for the map r ^ Hr- As [Hr + a] ha [Hr + a 


is contin¬ 


uous operation we see that Za{H^ G C'([0,T],U). It remains to show 
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that the norm estimate holds. Having e G (0,1), 
and a = 2e^ we have 


Mg - (i/ + a) 


H-H 


a 


C([0,T],V) 


< e, 


< Mg-. 

c([o,T],y) 2 


c([o,r],y)' 


= 1 and Za{H^ is the map 


r I—>• { Hr + a 


Thus Tf( Mg - {H + a) 


Let r G [0,T]. We denote 

Ha^r = {Hr + «), H^^r = {Hr + «), E = Ha,r “ H^^r- 
As Hr is positive semidehnite we have 


(30) 

Moreover 


\H. 


-i| 


y <a 


-1 


Thus 

(31) 


H-^ - H-^r = {[I + H-i-E] - /)M-i 

\\{Ha,r)-^E\\ 


< 


ly 


y- l-||(M„,,)-iE|| 


y 




We have ^ H Using (30) and (31) we have 


{H^^rr - {Ha,r) 


-1 


< 


y “ 1 - - 
2 


1 ||(^a,r) ||y < 2^2 — 2 e □ 


Let and B be as defined in ([^. We define 
(32) 5:C([0,T],F)^C([0,T]), 

S{Zr^){r) = {Z„(r)P^Pl,Pl)i2(o,2T)- 
Proposition 3. We have ||>S'||(^([o t]) — W’ 

Proof. As the maps r ha Pj.P 1 and r ha Za{r) are continuous, we have 
that S(^Za) G (^([OjT]). Let r G [0,T]. We have 


IP II <1 IIPIII^ = _ 


and therefore 




< — 

' - 3 

Zo.{r) 


□ 


y 
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Lemma 1. Let c G V{A). There is C > such that for all r > 0 
and p G H^{M) satisfying supp{p) C M(r) there is f E Sr such that 
u^{x,T) = p{x) and 


(33) ll/ll 

L2(o,2T) — C\\p\\ 

We recall that M(r) is defined by (10) and Sr is defined by 0 - 


Proof. Let us consider the wave equation with time and space having 
the exchanged roles 

(34) {dl - c{x)-‘^d^)u{x, t) = 0, (x, t) G (0, x(T)) x (0, T), 

u{x,T)=p{x), a:G[0,x(T)], 

u{x{T)S) = dxu{x{T),t) = 0, t G (0,T). 


By |37] the solution of ([34|) satisfies 


(36) 


■ )II_H'1(0,T) ^ C\\v\\ H^{M{T)) ■ 


If suppip) C M(r) then sMpp(n(0, ■)) C [T — r, T] and ^(a;, 0) = 
9tn(x, 0) = 0, when x G [0,x(L")], by finite speed of propagation. We 
choose fit) = M(0,f). □ 


Let fa,r be as in (12) and define 
(36) Sa G C([0,T]), Sa{r) := {fa,r, B1)]^2 ^q 2T)- 


Lemma 2. Let a G (0, min(l, Let V be as defined in (IS). 

Then there is C > 0, independent a, such that 

II-Sq — 14ilc'([o,r]) — Ca*. 

Proof. Let r G [0,T] and 5 > 0. Let us define ws G H^{M) 

( 1, if a;G(0,x(r)), 

W 5 (a;) = < jf X E[x{r),x{r) + d], 

[ 0, if a; G (x(r)+5, cx)). 

Using c(a;) > Cq we have 


(37) 


- 3 (^ 2 ' 


When 6 G (0, min(l, ^^^)) we have 


(38) 


wsfni^M) < X{T) + 


5 1 3 
3 5 - 5 
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Below C > 0 denotes a constant that may grow between inequalities, 
and that depends only on m, Co, Ci, L. Lemmagives us fs for which 
u^^{x,T) = ws{x). Thus (38) implies 

C 


II/5|Il2(0,2T) ^ . 


(39) 

Let f ^ Sr- We define 

Cra^rif) T O II / IIL2 (QgT) 


Using ( |37[ ) and ( |39| ) we have 

II 112 2 

(41) Goi^rifs) ||^'5 IIl2 (M W) ^ 11'^'^ 


(5 C 

L2(o,2T) — Q 


Functional ( [40| ) and the functional defined in Theorem]^ have the same 
minimizer fa^r- Using ([^, (21), and ([3^ we have 


l®« ~ ^llc(fori) = sup \{fa,r,Bl)L 2 ^^o^ 2 T]) -V{r)\ . 

re[0,T] 

= sup |(m-^“’’"(T), l)i2(Af;dV) — (Ijup), l)L2(M;dU)| 

re[0,T] 

<C sup \\uG-{T) - lM{r)\\l2(^M-,dV) - ^ <^«t(/«t) 


re[o,r] 

< C sup Ga,r{fs) 

re[0,T] 


re[0,T] 


Using (41) and choosing 6 = we have 

|2 


ka ^llc([o,r]) — Ga^. 


□ 


Lemma 3. Let V be as defined in (18). Let v be as defined in (19). If 
c G HiA) then there exists m > 0 s.t. 

ll^llc2([o,r]) < ^ ll^llc3([o,T]) < 


Proof. Equations (17), (18), (19), and (20) with the chain rule and the 
formula for the derivatives of inverse functions give us the result. □ 


For small h > 0 we consider the partition 

(0, T) = (0, h) U [h, 2h) U [2h, 3h) U ... U [Nh — h, Nh) U [A^h, T), 

where iV G N satisfies T — h < Nh < T. We define a discretized and 
regularized approximation of the derivative operator dr by 

(42) Z1;,:C'([0,T])^L'-(0,T), 
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Sqjh) 

h 2 

_ ) Sa{jh+h)-Sa{jh) 


if re (0, h), 


Dn{s^){r)={ if re[ 3 Kjh + hl 

Sq{T)-Sq(Nh) ^ re[iV/l,T). 

Proposition 4. Let B > t) and e e (0, min(^, —r))- Let a = BL^, 

/33 P^xiT)'^ 

h = 62 , V be as defined in (18) and let be as defined in (36). Let us 
assume that'sa e Sc([o,T])(sa, e). Then 


\\Dh{Sa) 

where C is independent of a and 5^. 


< CeT 


Proof. Let r e [jh,jh + h). Using (42) we have 

Sa{jh + h) - Sa{jh) 


Dh{sa)ir) - drVir) 

< 


— drV{r) 


Sa{,jh + h) - So,{ih + h) ^ Sc,{jh) - Saijh) 


h 


h 


+ 


Sa{jh + h) -V{jh + h) ^ V{jh) - s^Uh) 


h 

vUh + h)-vUh) 

h 


— drV (r) 


Lemma 


1^ gives us ||U||(^3qq < m. When r e [jh, (jh + h) there is 

^ e (jh, jh + h) such that 

U(jh + h)-U(jh) 


(43) 


h 


drV{r) = drV{0-drV{r) 


< hm. 


Using (43) and Lemmawith assumption we get 

2e 2Ca3 

D,{I^){r)-drV{r) 

h h 

Let us choose h = and a = Then 


+ hm. 


(44) 


Dh{sa){r) - drV{r] 


< Ce2. 


The proof is almost identical when r e (0, h) or r e [iVh,T). Note that 
the right hand side of (44) is independent of r. 

□ 

Let Cq and Ci be as in ([^. Let G L°°(0,T) and we dehne 

d-. if h„(r) < 
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We define 

(45) 

W : L“(0,T) ^ L~(M), W(fc„)(r) = ( ^ ^ (0,^), 

^ ^ ' V if r e [T, oo). 


Proposition 5. Let V be as defined in |ig| ) and v be as defined in 
(19). Let us assume that ka € BLrxi^QT){drV,e). Then 


W{k^) - V 


L°°{M) 


< C^e. 


Proof. For all x G M, we have 0 < Cq < c(x) < Ci. Let r G (0,T) 
and let assume that < ka{r) < Cf^. Using (19) and (20) we have 
0 < ^ < drV{r) < Then 


(46) 


1 

1 


ka{r) — drV{r) 

ka{r) 

drV(r) 


ka{r)drV (r) 


< Cle. 


In the case when r G (0, T) and k^ir) < Uf ^ or ka{r) > we obtain 
similar estimates. Note that the right hand side of (46) is independent 
of r. When r >T the left hand side is identically zero. □ 


Let Wa G L°° (M) and we dehne 

(47) 


Co, if Wa{r) < Co, 


T‘^:L“(M)^L“(M), T‘^(iF„)(r) := <^ w„(r), if Uq < ih„(r) < Ui, 


Cl, if Wa{r) > Cl. 


Let Wa G L°^{M) and we dehne 

(48) 


T : L~(M) ^ C{M), T{wa){t) = / Waiffidt'. 

Jo 


Using (47) and (48) we see that T o T'*’(w;q,) : M —)■ M is bijective as a 
function of t. Let us denote x = T o'^'^[wa) and x~^ = 

We dehne the sixth operator by 

(49) 

_ j Wa ox if X E [0,L), 

1, if X G [L, oo). 


$ : L“(M) ^ L“(M), <l>(w„) = 


Proposition 6. Let e > 0 and v be as defined in (19). Let us assume 
that Wa E B 1^00 (M){v,e). Then 


|<h(w;„) - cl 


L°°(M) 


< Ce. 
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Proof. Let us denote t = x and t = x Let x G [0, L). Using 
(20) and (49) we have 

\^{wa){x) - c(a;)| = \wa{t) - v{t)\ < \wc,{t) - v{t)\ + \v(t) - v{t)\. 
Lemma 1^ gives us ||u||c'2(ot) ^ and we have 

(50) \v{t) — v{t)\ < m\t — t\. 

Using ([^ and (20) we see that 0 < Co < v{t) < C\ and we have 

(51) Co|f-t| < I [ v{t')dt'\ = \x{t) -x{t)\- 


Having x(t) = x = x(t) and using (20) and (48) we see that 

(52) 

\x{t) - xit)] = \x{7} - x{7)\ = I [ (n(t') - Wa{t'))dt'\ < x~\L)e. 


Using (50),(^, and ( [5^ we have 

(53) 


l>i>(fi;„)(x) - cwi < (i + 


e. 


Note that the right hand side in (53) does not depend on x. When 
X G [T, cx)) the left hand side in identically zero. □ 

Proof of Theorem Let 

/ N r 1 1 

(54) Co = mmjl, —, 


2T’2xT9’2¥T9x(T)i 


}• 


Suppose that A G By{A, e) We denote H = HA and H = HA. Using 
Proposition we get 


H - H 


< 2Te. 


C{[0,T],Y) 

We denote Zo = Za{H) and Z„ = Za{Hy We have AT G i3c:-([o,T],y)(-ff, 2Te) 
and e G (0,min(l, ^)). Proposition [2 with P = | gives us 


Zr, — Zr, 


c([o,r],y) 


< =■ €i 


since a = 2 P+^TPeP = 2 ^Tiei. 


We denote Sa = SZa and Sa = SZa. We have G i3c([o,r],y)(ZQ,, ei). 
Proposition gives us 

2-2p 
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We denote ka = Dh{sa)- We have Sa € ^ci[o,T]){.Sa,e 2 ) and €2 
fo,niin(^, —r)V Proposition 0 with B = gives us 


- drV 


L°°(0,T) 


2~P 1 

< =: 
35 


es, 


e 


since a = ^ 2 ^Tiei. 

We denote Wa = W{ka). We have ka € i3i;,oo(o,T)(^rh", 63 ). Proposition 
1 ^ gives us 


\\Wa - v\\lo.(m) < Clc‘^T‘^ ^62 P =: 64 . 

32 

We denote = <h(u;a). We have Wa G 64 ). Proposition]^ 

gives us 


Cry C 


L^{M) 


< CCl‘^T^-Pe'^-P 
32 


Let e G (0,eo). Using (22),(28),(32),(42),(45), and (49) we dehne 
(55) 'kla{e) : U —)■ X, 

'Tlaie) = ^oWoDhoSoZaOH, 
and we have an estimate 


^a( 6 )(A) - c 


< Ceis. 


X 


□ 


Appendix A: The direct problem 


For the convenience of reader, we give the proof of the following, 
quite well known result, for continuity of the direct problem in our 
setting. 


Theorem 4. Let c G 'D[A) and f G L^(0,2T). Then the boundary 
value problem 0 has a unique solution G iL^((0,2T) x M). The 
operators A and Ut, defined in 0 and are linear and bounded 

operators, and the direct map A, defined in is continuous. 

Proof. Let consider the wave equation ([^. When c = 1 on M we 
denote the solution by Uq and have 


ul{t,x) = h{t — x), h{s) 


Jo f{t)dt, t > 0, 
0, f < 0. 


Notice that / h-)■ Mq is continuous from L^(0,2T) to iL^((0,2T) x M). 
Let us now show that the same is true for / 1 —>■ U'f. 
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We choose 'ijj G C°°{M) such that -0 = 1 near a: = 0 and c = 1 in the 
support of ijj. The commutator A = [8“^, ip] is a first order differential 
operator, whence Au^ G L^((0,2T) x M) for / G L^(0,2T). Let w be 
the solution of 

{d^ — c{xYd‘^)w{t,x) = Aul in (0, 2T) x M, 
d:,w(t, 0) = 0, 
w\t=o = dtw\t=o = 0. 


Then w G if^((0,2T) x M), see e.g. and U'f := ipu^ — w G 

H^{{0,2T) X M) is the solution of Indeed, c^dlip = ipd^ + A, 

where ip is interpreted as a multiplication operator. Here we are using 
the fact that c = 1 in the support of ip. Thus 

{dt - c^dDu^ = ip{d^ - dl)ul + Aul - Au^ = 0 . 


As "0 = 1 near x = 0, we see that u satisfies also the boundary condi¬ 
tions in ([^. 

Let us now suppose that / G C'“(0,T). Let be solution for the 
boundary value problem in ([^. Let c be as defined in Q. Let x G M 
and we define 

(56) k G k{x) = exp{ | ) 

^CyOC ] 

and 

(57) 

G : C‘^{M X (0,2T)) ^ C{M x (0,2T)), 

G{u) = k~^(^l — c^dljku = — c^dl — 2c^k~^dxkdx — c^k~^dlk^u. 

Let X G M and define 

(58) 0(x) = f c{x')~^dx'. 

Jo 

Let us denote x = (p{x) and define 


(59) 


U^{x, t) = U-^{(p{x),t) 


u^{x, t) 
k{x) 


Using (56), (57), (58), and (59) we see that {xA) is a solution of the 
boundary value problem 


(60) [8“^ — d~ + q{x))u^{xA) = (x, f) G M X (0, 2T), 

8iU^{QA) = 8 ^u^{LA) = Q, fG(0,2T), 

m7(x, 0) = 8tU^ {x, 0) = 0, X G M, 
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where 

(61) 


q(x) = -c (iji (x))k (4> (x))d^k(ili (x)) 


We define Xqf = u\x=q- Let us consider two velocity functions ci and 
C 2 , and let qi and q 2 be the potentials corresponding to ci and C 2 via 
formula ( [6l| ) . Using (59) we have 

(62) ||Aci ~ ^C2 il£(L2(0,2T)) — 11^91 ~ ^92 il£(L2(o,2r)) • 

Let us denote by and the two solutions with respect to potentials 
qi and q 2 for the problem ( [^ . Let us define w(x,t) = — 

Then w is the solution of 

(63) {d'f-dl + qi{x))w{x,t) = F{x,t), {x,t) e M x {0,2T), 

d^w{0,t) = 0, d^w{L,t) = 0, fe(0,2T), 

w(x, 0) = dtw(x, 0) = 0, X E M, 

where 

(64) F(x,t) = {qi(x) - q 2 (x))Pg^(x,t). 

By |SZ] we have an estimate 

(6b) ll'^llRi(Mx(0,2r)) — ^ ll-^llL2(Mx(0,2r)) • 


Let T = 0(L) and let supp{u^J C [0,T]. Using (64) we have 

(66) I|-^IIl2(mx( 0,2T)) ^ Ikl “ ^2||2,cx>(0,f) Ih^lUpMxCOgT)) ' 


Because / i-A is continuous from L^(0,2T) to H^{M x (0,2T)) we 
have 


(67) 


\u 


f 11 

<}2 llRi(Mx(0,2r)) 


< C 


L2(o,2T) • 


Using (|^ we see that 


( 68 ) 


WKf -KfWmopT) 


u 


f 11 

92 II_H- 1 (Mx( 0 , 2 T)) ■ 


Using ( |66[ ), ( |65[ ), ( |67[ ), and ( |68[ ) we have 

(69) 


Using 
and de: 


11^91 ^92 ll£(L2(0,2r)) — ll*?! '?2ilioo(o,T) • 

we have q{x) = —c^{x)k~^{x)d‘^k{x)\x=(j,-i{x)- Let x E M 


hi{x) = —ci{x)k^ {x)d^ki{x), h 2 {x) = —C 2 {x)k 2 {x)d^k 2 {x). 
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We have 


\hi{x) — h 2 {x)\ = \cl{x)k 2 ^{x)dlk 2 {x) — cl{x)k^ ^{x)dlki{x)\ 
< \c\{x) - cl{x)\\k];^{x)\\dlki{x)\ 

+ - k^^{x)\\cl{x)\\dlki{x)\ 

+ \dlki{x) - dlk2{x)\\cl{x)\\k^^{x)\ 


Using (56) we can bound each of these tree terms and have 
(70) ||gi — q 2 \\ x,oo(o,T) ^ ||ci — C 2 |Ic 2 (o,l) • 

Now (62), (68 ), and (iTOj) imply that A is continuous. 


□ 
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